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A detailed study of the various cosmological aspects in massive gravity theory has been 
presented in the present work. For the homogeneous and isotropic FLRW model, the 
deceleration parameter has been evaluated, and, it has been examined whether there is any 
transition from deceleration to acceleration in recent past, or not. With the proper choice 
of the free parameters, it has been shown that the massive gravity theory is equivalent to 
Einstein gravity with a modified Newtonian gravitational constant together with a negative 
cosmological constant. Also, in this context, it has been examined whether the emergent 
scenario is possible, or not, in massive gravity theory. Finally, we have done a cosmographic 
analysis in massive gravity theory. 
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I. INTRODUCTION 


In classical field theory, a basic question is whether it is possible to introduce a consistent 
extension of general relativity for massive graviton, or not. In fact, a small graviton mass may 
have significant physical interest, notably, for the cosmological constant problem. An attempt 


a. 


They 


for a covariant theory was initiated long back in 1939 by Fierz and Pauli (FP) 
introduced a quadratic mass term m 2 (h fiu h^ l/ — h 2 ) for linear gravitational perturbation 
to the action, and, as a result, there is a violation of gauge invariance in general relativity. 
Also, the linear theory with FP-mass does not recover General Relativity in the massless limit 
m —> 0, a contradiction with solar system test due to the van Dam-Veltman-Zakharov (VDVZ) 
discontinuity ( 2 ,( 3 ]. This discontinuity can be overcome by introducing nonlinear interactions 
with the help of Vainshtein mechanism [4]. But elimination of VDVZ potential led inevitably to 
the existence of Boulware-Deser (BD)-ghost [Jj], and, as a result, the theory becomes unstable 

Jfl- 

Then after a long gap, very recently, a two parameter family of nonlinear generalization of the 
FP theory was formulated by de Rham, Gabadadze and Tolley (dRGT) [10, ll[. This model of 
massive gravity is in the context of extra dimensions, and it modifies general relativity at the 
cosmological scale. Here, the BD ghosts are eliminated in t 
perturbation by constructing a covariant nonlinear action 


12 


re decoupling limit to all orders in 
-15 ] (see ref. [l(| for a review), and 


is free from tames issues which have plagued earlier attempts at a given mass to the graviton. 
In the context of studying the nature of perturbations on the given backgrounds, Fasiello and 
Tolley [ll] have investigated the Higuchi bound for classical stability of dRGT theory with the 


FLRW reference metric, and, in contrast to earlier investigations 


18-21], it is found that the 


bound is independent of the precise form of matter. Also, it has been shown [17| that the free 
parameters do not change the qualitative picture of the Higuchi-Vainshtein analysis except for 
a specific cosmological epoch. Furthermore, it is found that the Higuchi-Vainshtein tension can 
not be relaxed in the parameters space for stability of the dRGT theory if both the reference 
and the dynamical metrics are FLRW, i.e., the cosmological solutions (in massive gravity) with 
a general FLRW reference metric are not phenomenologically viable due to the tension with an 
operative Vainshtein mechanism [pil ]. 

At present, there are wide classes of applications of dRGT theory, particularly, in cosmology 


22-[32|], in Black holes, and, in spherically symmetric solutions 


bimetric gravity 


33-1401] and its relation with 


41-149]. In the present work, various cosmological aspects have been studied in 
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the perspective of the massive gravity theory. The paper is organized as follows: Section II deals 
with basic equations in massive gravity theory. Deceleration parameter has been evaluated, and 
possible transitions have been studied in section III. Section IV discusses on the equivalence of 
massive gravity and Einstein gravity. In Section V, we proposed an emergent scenario. Section 
VI is devoted for a detail cosmographical analysis. Finally, in Section VII, we have briefly 
discussed our results. 


II. BASIC EQUATIONS FOR MASSIVE GRAVITY 


Massive gravity theory has an effective field theoretic prescription given by Einstein gravity 
together with the covariant FP mass term. In dRGT model, the action is a functional of the 
metric g^vix) and four spurious scalar fields 4> a (x), a = 0,1,2, 3, which are introduced to give a 


manifestly diffeomorphism invariant description 


]. One defines a covariant tensor H^ as: 


9[iv — dv(t> Tjab T H[ivi 


(1) 


where g a b = diag(— 1, +1, +1, +1), the usual Minkowskian metric is defined in the co-ordinate 
system defined by cf> a, s. So, gravity in this formulation is described by the tensor H gu propagating 
on Minkowski’s space. In the unitary gauge, all the four scalars 4> a (x) are frozen, and, equal to 
the corresponding space-time co-ordinates cf) a (x) = However, sometimes, it may be helpful 


to use a non-unitary gauge in which </> a (x)’s are al 
density for massive gravity can be written as [10, 


owed to fluctuate. So, a covariant Lagrangian 

11 ] 


L = [R + m 2 g U(g , H)} , (2) 

where R is the usual Ricci scalar, m g is the mass of graviton and U, the nonlinear higher 
derivative term for massive gravity is a potential for the graviton which modifies the gravitational 
sector. The explicit form for U is given by 


U — U2 T CV3D3 + 


( 3 ) 
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where, 0 : 3 , aq are dimensionless parameters, and, Uf s {i = 2,3,4) have the expressions 


U 2 = 

[k] 2 -[k% 

(4) 

U 3 = 

[k ] 3 — 3[k][k 2 ] + 2[k] 2 , 

(5) 

U A = 

[k] 4 — 6[kt][k] 2 + 8[k 3 ][«;] — 6[k 4 ]. 

(6) 


Here, is defined by the four Stueckelberg fields cf> a (to restore general covariance [0j]) as 


= K ~ V9^ fabdA a du(!> b , (7) 

and the square bracket stands for the trace of the corresponding quantity within it, i.e., 


[*] = [« 2 ] = »£«£> 


and, so on.. Also, we write 


( 8 ) 


= d^d^fab, (9) 

as the symmetric (0, 2) tensor defined by the Stueckelberg fields (f> a . Normally, the reference 
metric f a b is chosen as the usual Minkowski metric g a b- Now, for convenience, if we choose the 
unitary gauge cj) a (x ) = then the metric tensor g /JV stands for the observable describing the 

5 degrees of freedom of the massive graviton. 

III. MASSIVE GRAVITY IN FLRW MODEL 

The line element for Friedmann-Lemaitre-Robertson-Walker (FLRW) geometry of arbitrary 
spatial curvature is given by 


ds 2 = g a pdx a dx^ = — N 2 (t)dt 2 + a 2 {t)^/ijdx l dx\ 
with spatial metric 7 ^ in terms of spherical co-ordinates as 


7 i jdx l dx J 


dr 2 

1 -Kr 2 


( 10 ) 


+ r 2 (d0 2 + sin 2 ddcf) 2 ), 


( 11 ) 
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where K = 0, +1, or —1 refers to flat, closed, or, open models of the universe respectively. 


In the present work, we choose the reference metric f ab as the de Sitter metric. The reason 
behind this choice in contrast to the usual Minkowski metric is that one can easily be able to 
formulate the flat, open, or, closed cosmologies by suitable slicing of the de Sitter model. Also, 


this choice of reference metric eliminates the problem of ‘no-go’ theorem 
model can be written as 


22|. The metric in this 


f ab dfd4 = -dT 2 + b 2 K (T) lij (X)dX i dX^ 

(12) 

b 0 (T) = e HcT , 

(13) 

b-i(T) = H~ l smh(H c T), 

(14) 

b\(T) = H~ l cosh(I7 C T). 

(15) 


It should be mentioned that in the limit, H c —> 0, one recovers Minkowski’s metric for flat 
and open cases: &o = 1,6 —i = T, while the later case corresponds to the Milne metric for flat 
geometry. 

Further, specifying the Stueckelberg field as 


cjP = T = f{t ), </>* = X' = x\ 


(16) 


one sees that the cosmological symmetries are satisfied and from Eq. ©, becomes a 
homogeneous and isotropic tensor of the form 


= dia S (~f 2 ’ b K (/(*))7*7) • 

So, from Eq. 0, the elements of k matrix have the simple form as 


(17) 


«o = l-C/^, ( 18 ) 

«$ = (! (19) 

J a J 

4 = 0 , ( 20 ) 

4 = 0 , ( 21 ) 
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where denotes the sign of ‘/’. Thus, in the Lagrangian, the non-linear higher derivative 
term for massive gravity (denoting the potential for graviton) becomes 


Lmg — V—gU ( g , H) — y/—g(U 2 + ot 2 ,U% + (X 4 U 4 ), (22) 


or, equivalently, 


Lmg — (a — bx{f)) [A r {a 2 ( 4 a 3 + 0:4 + 6 ) — a( 5 a 3 + 2«4 + 3 )bx{f) + («3 + a 4 )bj<(f)} 

— C//{ fl2 (3 + 3a3 + 0 : 4 ) — a(3a3 + 2a±)bK(f) + a 4^(/)} • (23) 


So, the variation of this Lagrangian with respect to ‘/’ gives the equation of motion for f(t) 


as 


[(3 + 3«3 + « 4 )a 2 — 2(1 + 2a 3 + 04 ) 067 ^ (/) + (03 + aq)^(/)] — C fb^if)] — 0- (24) 

From the square bracket, we have 


M/W) = /i±o(t), 


h± = 


1 + 2«3 + Ct 4 zt \J 1 + 0:3 + Q ! 2 — OL 4 


CV3 + 0:4 


(25) 

(26) 


So, one can obtain f(t) from Eq. (125]) . provided, b^ is invertible. It should be noted that, in 
case of Minkowski’s reference metric (i.e., f a b = 770 , 5 ), we have &o(/) = 1 for flat case, and hence, 
there is no longer any solution, while, for open case (where b-i(f) = /), there are two branches 
of solutions. On the other hand, from the remaining part of Eq. (1241) . we have 


CM) = (27) 

and as before, non-trivial solutions are possible, only if ‘ 6 ^’ is an invertible function. In this 
case, there is no analog of Minkowskian reference metric (as there exists no solution for both flat 
and open cases). However, with the choice of bx from Eqns. (usd, m, and (fl5l) . it is possible 
to have an explicit solution for /(f). In particular, for flat case (i.e., K = 0) (also assuming, 
/ > 0 ) we have 
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f(t) = H c l ln 


rr = 1 a 

JVa’ 


^ iL(f)a(f) ^ 


(28) 

(29) 


where H is the Hubble rate of the FLRW metric, and H c that of the reference metric (here, 
de Sitter metric). 


IV. FLRW COSMOLOGY IN MASSIVE GRAVITY THEORY AND THE POSSIBLE 
TRANSITIONS OF THE DECELERATION PARAMETER 


For the present FLRW metric, the usual Einstein-Hilbert (EH) term takes the form 


L EH = -^ + 3KNa, (30) 

together with which we have an arbitrary matter Lagrangian L m that describes ordinary 
cosmological matter. Thus, the variation of the total Lagrangian with respect to the lapse 
function N (which sets to unity in the following) gives first Friedmann’s equation 


3 ^ H 2 + — 2 ^ — 87T G(p m + p g ), (31) 

where p m denotes the ordinary matter energy density, and, p g is considered as an effective 
energy density due to massive gravity action as follows 

mP 1 

Pg = 8 jT Q a 3 ( b x{f) ~ a) [(6 + 4a 3 + a 4 )a 2 - (3 + 5a 3 + 2 a 4 )ab K (f) + (a 3 + a 4 )b 2 K {f)\ . (32) 

Further, the variation of the total Lagrangian with respect to a(t ) yields second Friedmann’s 
equation in the form 


2 H + 3 H 2 H—g — —87r G(p m +p g ), 
a z 


(33) 


where p m is the thermodynamic pressure of the cosmic matter and the effective pressure p g 


takes the form 





Pg 


m n 


{6 + 4cH3 -(- Q 4 — (3 -(- 3o?3 CX 4 ) f 


8 irGa 2 . 

+{1 + 2 0^3 + 014 — ( 0:3 + cxa) 
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2{3 + 3a s + ct 4 — (1 + 2a 3 + 0:4 )/}ofrftr(/) 

(34) 


As we are restricted to flat FLRW model, so, using the solution (|28p for fit), the expressions 
for p g and p g are simplified to 


Pg 


Pg 


m g 

8 v rG 


~Pg + 


H H 2 H ^ 

- 6a + 9/3 3^ +31^ 


m„ 


H 


87 rG HH r 


H 


-H 2 \ 


-3/3 + 27 — - 34-. 




if 2 / 


(35) 

(36) 


Note that, the present massive gravity theory contains three free parameters m g , a 3 , and 0 : 4 ; 
and the parameters a,/ 3,7 and 4 in Eqns. (13511 and (1361) are not independent, rather, they are 
related to 03 and 0:4 by the following relations: 


ol — 1 + 2 a 3 + 2 a 4 , 

(37) 

/3 — 1 + 3^3 + 4q(4, 

(38) 

7 = 1 + 6 a 3 + 12a 4 , 

(39) 

5 = a 3 + 4a 4 . 

(40) 


Interestingly, if H(z ) = H c , then there is no contribution of massive graviton to the energy 
density (i.e., p g = 0 ). Using the density parameters, Friedmann equation (I3TT) can be written as 
(with K = 0) 


n m + — 1, (41) 

where Q g =p g /p c ( p c = 3H 2 /8itG) is termed as the density parameter corresponding to the 
induced matter held in the massive gravity theory. Also, from the above Friedmann Eqns. m 
and (13311 . the expression for the deceleration parameter looks 


Q = 


\ (1 + 3Q m uJm) ~ + 9 f3-^ - 3q^ + 34-j|jj - ^mr c (“ 3 /3 + 2 7 ir c “ 3( +) 


h 3 \ 




+ nun ( _3 /3 + 2 7j 37 - 35 fj) 


(42) 
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where ui m = p m /Pm, is the equation of state parameter for the given fluid. 

FIG. 1 shows the graphical representation of q over the Hubble parameter (H) for the men¬ 
tioned set of parameters in the graph. We also mention that, for two very very small values of 
the graviton mass m g = 0.0000001 (dashed line), and m g = 0.00001 (dotted line), the graphs 
overlap with one another, and even for m g = 0.01, we can not distinguish between two graphs. 
It indicates the complete evolution of our universe from inflation to late time acceleration, where 
we have used the recent value of the density parameter [5l(. The right and left tails of the 


curves indicate the inflationary and the late-time accelerating phase of the universe, while, the 
matter dominated era corresponds to the middle portion of the graph, where q > 0. More¬ 
over, from the graph, we see that, if the graviton mass {m g ) is very very small, then the model 
approaches ACDM, however, if m g increases, then the null energy condition is violated at late- 
tirnes, and this may be potentially problematic for the stability of the theory. As the reference 
metric is de Sitter, and H c corresponds to its Hubble rate, so, H c would define where inflation 
may occur. However, in the present graphs, the chosen value of H c (= 1) falls into the matter 
dominated era. 

V. EQUIVALENCE OF MASSIVE GRAVITY THEORY AND EINSTEIN GRAVITY: 

COSMOLOGICAL SOLUTIONS 


In the previous section, we have presented effective Friedmann equations in massive gravity 
theory containing three free parameters. Here, we start with restricting the free parameters 0:3, 
a 4, such that, /3 = 5 = 0. As a result, we have a .3 = — aq = |, thus, we have a = and 

7 = -i 

So, the Friedmann Eqns. (1311) and (1331) are now simplified to 


0 3 o 3 m? n „ 

3 H 2 = 8 irGp m ~ ~m 2 g + ^f# 2 , 

9 

771 

2 H = —8 tv G(p m + Pm) + H, 


which can again be rewritten as 


(43) 

(44) 


3 H 2 = 8 'KGmPm - A, 

— 8 'TTG m (p rn T Pm )j 


(45) 

(46) 
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— ■ m =0.0000001 •••mm =0.00001 m =0.2 

_ g _ g _ g 


FIG. 1. The figure shows the cosmic evolution from 
inflation to late time acceleration (see Eq. (I42f) l for 
the following parameterizations: 03 = —1; aq = —2, 
c o m = 0, Cl m = 0.3086 and H c = 1. 


with G m and A as 


G 


m — 


A = 



(47) 

(48) 


Note that, Eqns. (I45p and (14611 are nothing but usual Friedmann’s equations in Einstein 
gravity with a negative cosmological constant which depends on the graviton mass, and also, the 
Newtonian gravitational constant is modified by the free parameter m g . It should be mentioned 
that, if the graviton is assumed to be massless, then we get back usual Newton’s gravitational 
constant, and the cosmological constant vanishes. In the literature [52-56], it is claimed that 


negative cosmological constant has a significant role in describing the evolution of the universe. 
The low energy limit of supersymmetry prefers a negative cosmological constant implying an 
Anti de Sitter (AdS) cosmos [57|]. Also, the problem with negative cosmological constant in 


supersymmetry theory can be resolved using Einstein-Cartan theory 


58 ]. It was speculated by 
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Biswas et al. 


59( |. that our universe was begun with a negative cosmological constant. However, 


Kallosh and Linde 


60] claimed that in presence of a negative cosmological constant, our universe 


may recollapse in future. 

Now, we shall present some cosmological solutions in massive gravity theory considering its 
equivalence to Einstein gravity with negative cosmological constant. 


A. 


When matter in the form of a perfect fluid with constant equation of state 


If the matter is chosen as a perfect fluid with constant equation of state (EoS) as, p m = uj m p m -, 
Lo m is a constant, then the Hubble parameter and the scale factor can be obtained by solving 
the field equations (1451) and (1461) as 


H = 



|A| 


tan 



|A| 


(A - e 0 t) 


cos 


a = B 

Also, the deceleration parameter (q) takes the form 



>-«*)' ' 


(49) 

(50) 


Q = -i + <+ 


1 + cot 2 



(51) 


Here, A, B are the constants of integration, and, eo is the slow roll parameter in the case 


when A = 0, i.e., eo:= — ( -jp 


A—>-0 


= |(1 + u! m ). The above solutions represent a model of the 


cyclic universe compared to the power law expanding solution in standard cosmology provided 
(8 tt Gm ) _ A >Q, i.e., if the initial energy density of the fluid (= po) is large enough to 


compensate the negative cosmological constant 
Now, the expression of the Ricci scalar 


611 ]. and po depends on the graviton mass. 


R = 6 (H + 2H 2 ) = 2|A| 


(2 — eo) tan 2 


■L-^(A — e 0 t) ] - e 0 


(52) 


shows that R diverges at td = ^ 
key feature of the cyclic universe, i.e., a universe with successive expansion and contraction, 
and continues in this manner, and this l td depends on the massive gravity effect through the 
cosmological constant in Eq. (|48|> . Thus, this cyclic nature of the universe comes from the 
massive gravity effect. 


A + ./t?t(2 n + 1)5 , n € Z (set of integers), and it is the 
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B. When the fluid obeys the Equation of state (EoS) for Chaplygin gas 

The EoS for the Chaplygin gas is p = — —, where A > 0. Now, inserting this EoS in Eq. 


(1461) . and integrating we have 


H 2 = -8vr G m VA 


1 + 


D 2 
Aa 6 


- 1 + 


D 2 V 

AaP ) 


da 

a 


(53) 


Here, we consider the solution for large a. Therefore, using binomial expansion, Eq. (1531) 
gives a solution for the Hubble parameter as 


H = 


l47rD 2 G m \ 1 


3 A^ 


(54) 


and, consequently, solving Eq. (1541) . we have the scale factor as 


a 3 = 


\2'kD 2 G ti 

Ah 


t + E. 


(55) 


The deceleration parameter q takes the form 


6 


1 


Q ^ l + H 2 da) 1 + 1 3 V 8vr D 2 G m I a 3 ' 


(56) 


Here D and E in the above equations are the constants of integration. Hence, we have a 
power law expansion of the universe, and the model gradually approaches ACDM (q = —1) for 
large a. 

Further, for this perfect fluid model, the general held Eqns. (ED) and (1331) can be rewritten 
as (for simplicity, we assume, 8 nG = c = 1.) 


3H 2 (1 + 7^|) = p 0 a 3 ( 1+w -) + 3 m 2 g ^-2a + 3/3-^- + 5^—) , 

p 0 HH c (l+u m )a- 3 ( 1 +^ 


H = — 


2 HH C + m 2 g (-3/3 + 2 7 ^ - 3hg) 

where p m is eliminated by solving the matter conservation equation. Thus, evolution Eq. 
with the definition of the Hubble parameter 


(57) 

(58) 
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a = aH , (59) 

together form an autonomous system in the phase plane (a, H), provided, 1 + ui m / 0, i.e., 
the fluid model is not a ACDM model. Note that a and H are not constants rather they are 
constrained by the relation (1571) . The critical points of the dynamical system are given for the 
following two regions 

-k when the fluid is not in the phantom region, i.e., 1 + u m > 0. The critical points are 

i 

( \ 3(l+wm) 

(fi=0,a = «„); • < 6 °) 

* when there is phantom fluids, i.e., 1 + ui m = 0, 


{H = Hq, a = 0). 


(61) 


The parameter Hq in Eq. (1611) satisfies the cubic equation 


H 6 - 


{He + 7 m 2 q ) 


m? g 5 


H 2 + 


3 f3H 2 c tt 2 aHl 


-H — 


= 0 . 


(62) 


As Hq is the real root of the cubic Eq 
choice of the parameters a 3 ,a^,m g , and H c . 


5 5 

, so the number of critical points depends on the 


VI. EMERGENT SCENARIO 


The proposed cosmological scenarios to avoid the initial singularity (big bang) of the standard 
cosmology can be classified as bouncing universes, or, the emergent universes. In this section, 
we shall examine whether in the context of massive gravity theory, it is possible, or not, to 
have the emergent scenario which arises due to the search for the singularity free inflationary 
models in the framework of standard cosmology. Briefly, emergent universe is a model universe 
having no time-like singularity, ever existing and almost static behavior in the infinite past (t —> 
—oo). It is worthwhile to mention that, in 1988, Brandenberger and Vafa [62(] showed that in 
the context of superstring theory, a simple model can resolve this big bang singularity. After 
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that, Ellis and Maartens 


63], and Ellis et al. [64| discussed the same scenario of the universe 


to remove this big bang singularity. Subsequently, Mukherjee et al. 


65] found some solutions 


for Starobinsky model having features of an emergent universe. Also, Mukherjee et al. 


66 ] 


formulated a general framework for an emergent universe considering “GR+ exotic matter”. As 
a result, several works were done on the emergent universe [67-711]. Recently, in the context of 


non-equilibrium thermodynamics, gravitationally induced particle creations exhibit the emergent 
universe scenario 


72]- 


In the context of massive gravity theory, we start with the Friedmann Eqns. (1451) and (14611 . 
These equations can be interpreted as a single fluid in Einstein gravity as 


3H 2 = 8irGp eff , (63) 

2H =-8irG(p eff +p eff ), (64) 


where 


Peff 


Peff 


Pm 


A 

8nG 


G 2 
p™ 877 ? 

n _ 22l)' 

l 1 2H?) 


The energy conservation equation 


(65) 

( 66 ) 


Pm T 3H(p m Pm) — 0, 


(67) 


now becomes 


Peff + 3H{Peff + Peff) — fl¬ 
it we now assume the equation of state parameter of the effective fluid as 


( 68 ) 


Peff = - jjPeff, 

where e > 0 is some constant, but depends on the mass of the graviton. 
(1681) with the assumption in (1691) we have 


(69) 

Now, integrating 
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o ,, = ?°e 3et 
" e ff n 3 C ‘ 


Further, using (1691) in the Friedmann Eqns. (1631) and ()64[) . we have 


(70) 


2 H = -3 H 2 + 3 eH, (71) 

which on integration gives 


H = 


3e 


3 + e 


(72) 


where r = 4r (t — to), and to is the constant of integration. Now, integrating the above relation 
once more, one gets the scale factor for the emergent scenario as 


x 3 

a \ 2 

— = 1 + 3e r , with ao = a(to). 

aoj 


(73) 


From the above solution in Eq. (1731) . we see that, as t —> —oo, a —> ao, and, for t <C to, 
a ~ ao, while the scale factor a grows exponentially for t > to- In particular, the above cosmo¬ 
logical solution has the following asymptotic behavior: 


(I) a — > ao, H— > 0, as t— > — oo 

(II) a ~ ao, H ~ 0, for t<< to 

(III) a ~ e H o{t-t 0 ) ^ f or to. 


Hence, it is possible to have an emergent scenario in massive gravity theory for a perfect fluid 
with variable equation of state (given in Eq. (1691) 1 depending on the mass of the graviton. 


VII. COSMOGRAPHY OF MASSIVE GRAVITY 

In cosmography analysis, the universe is assumed to be homogeneous and isotropic on the 
largest scale and no specific dynamical theory is assumed a priori. The scale factor is expanded 
in Taylor series with respect to the cosmic time as 
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n (f.^ — ^ + Ho(t — to) ~ 7^QoHo(t — to)~ + —joH^(t — to) 3 + —soH^t — to) 4 + 

^loHo(t - t 0 ) 5 + ^m 0 i^o(i - t 0 f + 0(\t - t 0 | 7 ); (74) 


where to is the present time, and, the suffix ‘O’ indicates the value of the quantity at present. 
The above coefficients of different powers of ‘t’ are defined as |3| 


d 1 cPa . 1 d 3 a 1 d 4 a 1 d 5 a 1 d 6 a 

a ’ ^ aH 2 dt 2 ’ ^ aH 3 dt 3 ’ S aH 4 dt 4 ’ aH 5 dt 5 ’ m ai7 6 dt 6 ’ 

and are conventionally termed as Hubble, deceleration, jerk (j), snap (s), lerk (l) and m 
parameters respectively. These parameters can be used to find the distance-redshift relation 
and hence different distances in the universe. Further, the sign of q indicates acceleration (—ve), 
or, deceleration (+ ve), similarly, a change of the sign of j (in an expanding universe) signals 
that the acceleration starts increasing or decreasing. 

In this context, we have studied the cosmographic parameters for three different values of 
the graviton mass ( m g ). 

FIGs. 2-5 show the variations of the four cosmographic parameters j, s, l and m against the 
Hubble parameter ( H). In FIG. 2, as there are several transitions of j, so the acceleration is 
fluctuating throughout the evolution as reflected from FIG. 1. 

VIII. CONCLUSIONS AND OUTLOOK 

In massive gravity theory, we have studied several cosmological aspects. The background 
metric is taken as flat FLRW metric, and the matter is chosen as perfect fluid with constant 
equation of state. The deceleration parameter has been calculated, and FIG. 1 shows the transi¬ 
tion from deceleration to present late time acceleration as predicted by the recent observations. 
Moreover, from FIG. 1, it has been found that the model will remain stable as long as the 
graviton mass is very small. However, the model gradually becomes unstable as the mass of the 
graviton gradually increases. Among the three free parameters in the massive gravity theory, 
by choosing a 3 and 04 appropriately, it has been shown that the Friedmann’s equations in this 
theory are equivalent to those in Einstein gravity with a negative cosmological constant, and, 
the Newton’s gravitational constant is modified by the remaining free parameter (the mass of 
the graviton m g ). Thus, the choice of the free parameters 0:3 and a4 changes the qualitative 
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FIG. 2. This is a study of the cosmo- FIG. 3. A study of the cosmographic parame- 
graphic parameter j described in the pa- ter s described in the paper for three different 
per for three different graviton mass m g . graviton mass m g . Parameters are a 3 = —1, 
Parameters are a 3 = -1, a 4 = 0, H c = 1 a4 = q, H c = 1 and oj m = 0 and Q rn = 0.3086 
and w m = 0 and Q m = 0.3086 [5l| . 511_ 




FIG. 4. This is studied for the cosmo- FIG. 5. It is a study of the cosmographic pa- 
graphic parameter l described in the pa- rameter m described in the paper for three 
per for three different graviton mass m g . different graviton mass m g . Parameters are 
Parameters are a 3 = -1, a 4 = r Jl H c = 1 a 3 = -1, a 4 = 0, H c = 1 and uj m = 0 and 

n m = 0.3086 [51]. 


and ui m = 0 and VL m = 0.3086 I51H . 
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: feature of the massive gravity theory, and it is in difference with the work by Fasiello and Tolley 
l?j]. The interpretation of negative cosmological constant has been discussed. A possibility of 
emergent scenario has been investigated in massive gravity theory. It is found that for a specific 
variable equation of state (related to the mass of the graviton) with negative pressure, it is pos¬ 
sible to have this non-singular aspect at early epoch which was also obtained in the superstring 


theory 


621 ]. Lastly, the cosmographic analysis has been done for this present gravity theory. 


The cosmological solution for perfect fluid with constant equation of state shows a periodic 
nature of the scale factor which may be interpreted as the possibility of a cyclic universe in 
massive gravity theory. This feature is also supported by the divergence of the Ricci scalar at 
infinite number of points. Also, it should be mentioned that, this cyclic nature of the universe 
which comes from the massive gravity effect can be realized in the superstring theory as discussed 
by Brandenberger and Vafa in Ref. 62]]. For Chaplygin gas model, the solution shows a power 
law expansion of the universe, and the model as expected approaches the ACDM model. Finally, 
the negative cosmological constant can be incorporated into the matter sector by introducing 
effective density and pressure in the Einstein equations. 
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